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Abstract 

Two particle correlations are used to extract information about the character- 
^ istic size of the system in proton-proton and heavy ion collisions. The size of 

the system can be extracted from the Bose-Einstein quantum mechanical ef- 
^ feet for identical particles. However there are also long range correlations that 

C shift the baseline of the correlation function from the expected fiat behavior. 

K- 5 A possible source of these correlations is the conservation of energy and mo- 

t~-. mentum, especially for small systems, where the energy available for particle 

^ production is limited. A new technique, first used by the STAR collabo- 

r— I ration, of quantifying these long range correlations using energy-momentum 

Q conservation considerations is presented in this talk. Using Monte Carlo sim- 

p ^ ulations of proton-proton collisions at 900 GeV, it is shown that the baseline 

D of the two particle correlation function can be described using this technique. 

> 1 Introduction 

CN Femtoscopy is a powerful technique that allows us to "look" directly into the 

interaction region in particle collisions. It was first used in proton-antiproton 
^ collisions by Goldhaber,Goldhaber, Lee and Pais in 1969 [1] but it was de- 

^ veloped by Hanbury Brown and Twiss in 1956 to measure the angular size 

'T^ of stars [2] and is therefore also known as HBT. With this technique the size 

^ and lifetime of the emission region can be measured to study the shape and 

^ evolution of its components. The measured radii and their dependence on 

^ pair momentum and multiplicity provide information about the level of in- 

teraction of the created particles, and can be used in particular to determine 
the freezeout volume in heavy ion and proton-proton collisions. 

The HBT effect arises from the Bose-Einstein quantum mechanical effect 
for identical particles. However, any given pair of particles can be correlated 
through other effects, like resonance decays, jets and mini-jets, or phase-space 
constrains from energy momentum conservation. In this paper the latter are 
studied, and it is organized as follows: a brief introduction to the femtoscopic 
techniques is given in section [2} the energy momentum correlations are ex- 
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plained in section [3] and their application to correlation functions are shown 
in section m 



2 Femtoscopy 

When dealing with two identical particles in quantum mechanics it is neces- 
sary to describe them with a symmetrized two-particle wavefuction: 

^(pi,ri,p2,r2) ~ e^(Pi-'^i+^'2-"2) + e^^Pi-'^+f^-'i) (1) 

Wavefunction symmetrization leads to an enhanced probability of measuring 
pairs of particles with small momentum difference: 

l^^l ~ 1 + cos((pi - P2) ■ (ri - ra)) = 1 + cos(Q ■ Ar), (2) 

where Q = Pi — P2 is the momentum difference and Ar = ri — r2 is the 
relative emission point of the particles, see Fig. [TJ 



Pairs with small relative momen- 
tum Q will be observed with 
higher rates because identical 
bosonic particles are likely to be 
found in the same state. This 
gives rise to a peak at low mo- 
mentum difference, and its width 
is inversely proportional to the 
size of the emission region. In 
particular for a Gaussian source 
p(r) ~ exp"'' /^"^ one obtains: 
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Figure 1: HBT effect. 



Q,(Q) = 1 + Ae-«'«^ (3) 

where A is the correlation strength. The experimental correlation function is 
defined as the ratio of the two particle and single particle distributions: 

Eqn. |3] assumes the only correlation between particles is bosonic. The AL- 
ICE collaboration measured Bose-Einstein correlations of pions at 900 GeV 
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[1], and it was observed that the extracted radii depend on what type of 
basehne is used in the fitting procedure. In this paper we are concerned 
about studying a possible source of non-femtoscopic correlations that can 
affect the measured radii, namely Energy-Momentum Conservation Induced 
Correlations (EMCICs). 

3 Energy-momentum conservation induced cor- 
relations 

The phase-space available for particle production is constrained by energy 
momentum conservation, therefore creating a correlation among particles. 
It has been shown by Chajecki and Lisa [5J, that these correlations can be 
quantified by applying energy-momentum constraints on the single particle 
distributions. This results in an expression for the energy-momentum corre- 
lation of two particles: 

where is the total event multiplicity. When femtoscopic correlations are 
present the total correlation function can be written as: 

C{Q) = <^femto{Q) X CeMCIc{Q) (6) 

where $ jemto is the femtoscopic or Bose-Einstein correlation. This expression 
implies that if all the quantities in Eqn. [5] could be measured, the EMCICs 
would be completely characterized and it would be possible to remove their 
effect from the experimental correlation function. However, in experiment 
A^, (pr), {vDi (E) and {E'^) cannot be measured, simply because not every 
particle is detected. Parametrizing the unknown quantities yields an equation 
that can be used to fit experimental data, as has been done in STAR for 200 
GeV p + p collisions [6] and is being studied in ALICE for 900 GeV p + p 
collisions |7j: 



C{Q) = l-Mr{pTi ■ pT2}-M2-{p,, ■ p,,}-Ms-{E, ■ E2}+Mr{Ei + E2}-- 

(7) 

The notation {X} represents histograms of two-particle quantities that can 
be measured in experiment and are binned at the same time as the numerator 
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and denominator of the correlation function in Eqn |4} The Mi parameters 
are related to the quantities that are not directly measurable in the following 
way: 



Ml = — ^ M 



N{vl) N{pl) 

1 [E) 
" N{{E^) - (E)2) ^ N{{E^) - (E)2) 

After a fit has been performed on the data, the physical quantities can be 
calculated using an additional equation for the energy of a characteristic 
particle with mass m^. in the system : 

(i^') = kv%) + (P^) + ml (10) 
The total multiplicity of the event is then calculated as: 

AT = ^3-^ -2^1" -^2" (11) 

The physical quantities can be used as a consistency check to tell if the fit 
results are reasonable. 



4 Baseline study of two-pion correlation func- 
tions. 

The shape of the baseline is studied using Monte Carlo simulations of proton- 
proton collisions at 900 GeV, with a data set consisting of lOM events. The 
data analysis used to obtain all the histograms mentioned above, is done 
using the AliFemto code. A transverse momentum cut is applied to ensure 
correct identification of pions in a detector like ALICE, i.e. 0.1 GeV <pt < 
1.2 GeV. The analysis also includes binning in four pair momentum bins 
kx = 1/2|pti + PT2I, whose edges are 0.1,0.25,0.4,0.55,1.0 GeV. 

A fit to the Qinv correlation function with no kx binning can be seen in 
Fig. [2j The behavior at low Qmv is from two track effects, known as merging 
and splitting. At large Qmv the long rage correlation is very pronounced. 
The fitting function used is Eqn. [6} including the Gaussian component to 
account for the underlying event at low Qinv The fit adjusts very well to 
the data, however, it was observed that the Mj values are not unique and 
depend on the initial values given to the parameters. This means that there 
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Entries 
N 


134955 
1.484 ±0.008 




Physical quantities 

N-1.S68 -f. 0.689 
<E> = 2.337 +- 0.010 
<E=> = 10.15+- 3.207 
<P^> = 5.147 +.2.280 
<P;>" 4.991 +.2.255 


lambda 

R 

M1 

M2 

M3 

M4 


0.4856 ±0.01 06 
0.2399 ±0.0018 
0.2477 ±0.01 33 
0.1277 ±0.0608 
0.1 358 ±0.0005 
0.31 75 ±0.0007 
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Figure 2: EMCIC fit to the Qinv correlation function. 

are many local minima, which is expected from a one dimensional function 
with 7 parameters. The physical quantities obtained from the Mj, also shown 
in Fig. [2| indicate that the fit is in a physical region. The average energy 
< E >= 2.34 ± O.OlGeV is the most significant here given the small error 
bar, and it is a reasonable value. The total event multiplicity is much lower 
than expected and the average transverse momentum is quite high but both 
have a large error bar. 



A way to constrain the parameter space is to use the kx binning mentioned 
above and fit all bins simultaneously with independent N, A, Rinv but the 
same Mj. This is possible since the Mj's are related to physical quantities 
which in this case are the same for all bins. The results of this procedure 
are shown in Fig. [3j The fitting function adjusts very good in all four fc^ 
bins, and it was observed that the convergence of the Mj values was more 
consistent when modifying initial values. As before the average energy < 
E >= 2.22 ± 0.026*6^^ is the most significant quantity, due to its small error 
bars, and it is again a reasonable value indicating a good fit. 

5 Summary and future work 

A new technique to quantify and remove effects of energy-momentum conser- 
vation correlations from femtoscopic effects has been used here to study the 
baseline of the Qinv correlation function. It was shown that the convergence 
of fits to the data is dependent on initial parameters but it explains the long 
range correlations from first principles. Future work on this topic includes 
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Figure 3: EMCIC fit to the Qinv correlation function simultaneously in 4 kj- 
bins. 

using the same procedure to fit correlation functions in the longitudinal co- 
moving system LCMS and in spherical harmonics, where there should be 
more constrains on the parameters allowing the fit to converge easier. 
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